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\ We consider a generalization of the interior Schwarzschild solution 

that we match to the exterior one to build global C 1 models that can 
have arbitrary large mass, or density, with arbitrary size. This is pos- 
sible because of a new insight into the problem of localizing the center 
of symmetry of the models and the use of principal transformations 
to understand the structure of space. 



1 Introduction 

X 

We consider in this paper a family of spherically symmetric, static models 
with bounded sources. To start with we shall consider the following reduced 
form of the line-element: 

ds 2 = -A 2 dt 2 + ds 2 , (1) 

with 

ds 2 = B 2 dr 2 + BCr 2 dtt 2 , dVl 2 = d9 2 + sin 2 9d(p 2 , (2) 

A, B, and C being three functions of r. This form of the line-element is fully 
adapted to the assumption requiring the existence of a global time-like inte- 
grable Killing vector field, as well as to the assumption of spherical symmetry. 
It remains invariant under the adapted time coordinate transformation 

i=Kt + K 0l (3) 
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where K is an arbitrary positive constant which is usually chosen such that 

lim A = 1, (4) 

9 — >oo 

and K is a second arbitrary constant that can be completely ignored. It is 
also invariant under a radial coordinate transformation 

f = f(r), (5) 

which may serve different purposes: mathematical simplicity or desired in- 
terpretation. 

Among the coordinate conditions that can be used to choose r one finds: 

i) The historical condition 

A 2 B A C 2 r A = 1. (6) 

This is equivalent to the coordinate condition used by Schwarzschild to derive 
his exterior solution. 

ii) The curvature condition 

BC = 1. (7) 

It was first considered for mathematical simplicity by Droste |IJ, Hilbert 
and Weyl [[| and it is the most used. 

iii) The isotropic condition 

C = B. (8) 

iv) The harmonic condition 

where the prime means a derivative with respect to r. This condition is 
equivalent to requiring that the three coordinates 

x 1 = r sin 9 cos (p, x 2 = r sin^siny?, x 3 = rcos# (10) 

be harmonic functions in the space-time defined by (H), i.e., 

1 „ / ^- j, 



Ax k = -j=d % (^gVd jX k ) =0, g = det (g af) ) . (11) 

(Greek indices run from to 3 and Latin indices from 1 to 3.) This condi- 
tion which leaves open the choice of an arbitrary constant was used, with a 
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particular value of this constant, by Fock [|]] who meant to give a particular 
meaning to this particular radial coordinate r. 
v) The Gauss condition: 

B=l. (12) 

In recent years one of us (LI. B) has introduced the quo-harmonic condi- 
tion HI: 



Ur>C)'. (13) 



B 

2r 

This is equivalent to requiring that the three functions (|T0|) be harmonic in 
the 3-dimensional space with line-element (|2|), i.e., 

Ax k = -^d t (yfgg ij diX k ) =0, g = det ( 9ij ) , g* = g'K (14) 

This new coordinate condition has been thoroughly discussed in 0, 0. 
It is intimately connected with the concept of principal transformation of 
a 3-dimensional Riemannian metric J5|. Although not every such metric 
possesses a principal transform, many do have it and in particular those with 
spherical symmetry. In this particular principal transform of (EH is a 

new metric with line-element 

ds 2 = § 2 B 2 dr 2 + ^ 2 BCr 2 dn 2 , (15) 

where $ and \1/ are two functions of r such that 

i) the Riemann tensor of (|15D is zero, 

Rijki = 0, (16) 

i.e., the metric is flat, and 

ii) the transformation from (|2T) to (O) is harmonic: 



(rj* - ry 9 jk = o ; (17) 

where f * fc and r*- fc are respectively the second-kind Christoffel symbols of (0) 



and ([15]). Since both ( |i"6"D and (|T7| ) are tensor equations as long as only 
coordinates that are adapted to the main Killing vector are considered, both 
conditions are intrinsic to the static character of the models. 

To solve for $ and ^ two methods can be used. The first one consist in 



using directly Eqs. (|T6|) and ([[7]) where these two functions are the unknowns. 
It has the advantage of allowing to use any coordinate condition that one 
wishes. The second method is based on the remark that if we consider the 
problem solved and we use Cartesian coordinates x % of (|15|), i.e., such that 

r}* = o, (is) 
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then from ( |TTD it follows that in this system of coordinates we shall have 

ry jk = o- (19) 

These conditions tell us that x % is a system of quo-harmonic coordinates 
of (@), i.e., solutions of fll4j) . /,From this remark it follows that if we solve 
first these equations for the functions 

x 1 = R(r) sin 6 cos (p, x 2 = R(r) sin 6 sin (p, x 3 = R(r) cos 9, (20) 

that is to say, if we first solve the single equation for the unknown func- 
tion R(r) 

r 2 R" + r (2 - rB^B') R' - 2B 2 R = 0, (21) 



then to solve Eqs. ( |TED and ([T7|) we just have to equate the line-element of 
flat space written in polar coordinates R, 8, ip with ([15]): 

dR 2 + R 2 dQ 2 = <S> 2 B 2 dr 2 + V 2 BCr 2 dtt 2 , (22) 

whence it follows that 

R 1 R 

^b' ^WW (23) 

Principal transformations were introduced in || as a generalisation to 
3-dimensional Riemannian metrics of one of Gauss's theorems according to 
which any 2-dimensional metric can be mapped conformally, locally, into an 
Euclidean space. In this paper they will play also an important role in the 
interpretation of the models to be presented in Sect. ||. 

These models will be required to satisfy the following conditions: 

i) A value of the coordinate r exists, say 7*1, such that for r > n the 
space-time model is a solution of the vacuum field equations 

S a p = 0, (24) 

where S a p is the Einstein tensor. It is therefore the exterior Schwarzschild 
solution. 

ii) For r < ri the space-time model is an interior solution with a perfect- 
fluid source 

S aj3 = KT a/3 , k = 8n, T af3 = (p + p)u a up + pg a p, (25) 

where p is the energy density, p the pressure and u a the 4-velocity of the 
fluid, tangent to the main Killing field. The solution that we shall consider 
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is a rather straightforward, but physically innovating, generalization of the 
interior Schwarzschild solution. This generalization is crucial to the goal that 
we pursue in this paper. Namely the possibility of constructing global models 
with arbitrary mass or arbitrary density, and arbitrary size. 

hi) On the border r = r\ we shall require the continuity of the three 
functions A, B, C as well as the continuity of the three derivatives A', B', C 
thus completing the construction of models of class C 1 . This requirement 
greatly restricts the choice of a unique global coordinate condition, i.e., being 
the same on both sides of the border. 



2 The exterior Schwarzschild solution 



The exterior Schwarzschild solution has been described using the coordinate 
conditions i) to iv) mentioned in the preceding section. We list below the 
corresponding expressions of the coefficients A, B and C, and the intervals 
of the corresponding radial coordinate r on which the metric is static, 
i) Historical condition:^] 



.4 



2m 



B = A- l B!, VBC = R/i 



where R(r) is the function 

R(r) = [r 3 + (2m) 3 
ii) Curvature condition: 



1/3 



.4 



2m 



B = A~\ 



iii) Isotropic condition: 

1 — m/2r 

= / 5 

1 + m/2r 

iv) Particular harmonic condition 



B = C 



BC 



m\ 2 
1 + Tr) 



re]0,oo[, (26) 



(27) 



r G ]2m, oo[ . 



r G ]m/2, oof . 



(2f 



(29) 



.4 



<r — m 

r + m ' 



B = A~ l , 



BC= 1 



m 
r 



r G ]m, oo[ . (30) 



1 Schwarzschild's original solution has been reconsidered recently in Ref. ||. 
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v) Gauss condition: 




A=Wl-— , B = l, C = R, re]D, oo[, (31) 
where D is the arbitrary constant in the following definition of R. 



R = Jl-^ + mln (-l + L + LM-^L )+D. (32) 

The lowest value of r, say re, in the domain of staticity, different for each 
of the coordinate conditions considered so far, defines a border, say B, which 
was used to be called the Schwarzschild singularity and that now is more 
often called the horizon of the exterior Schwarzschild solution. 

In the case of the original Schwarzschild coordinate, but also when us- 
ing the Gaussian condition and setting D = 0, one has tb = 0, strongly 
suggesting that Schwarzschild had succeeded, as it was his intention, in iden- 
tifying the most extreme source of his solution as being a point. This point 
of view it is still accepted in |1(| where the author claims that the non zero 



values of rs, suggesting that B is in fact a 2-dimensional surface, come from 
an inadequate choice of the coordinate condition. We want to comment that 
whether B is a point or a surface depends intrinsically on the geometry which 
is describing the structure of the space in the static reference frame. If the 
metric of this space is that with line-element (0), i.e., if one accepts that the 
physical distance d(xi,x 2 ) between two points of space is 

d (xi, x 2 ) = / ds, (33) 

J Xl 

then B is the border of the metric completion of the domain of staticity, 
and this border is unquestionably a 2-dimensional sphere. No coordinate 
condition can change that. On the other hand there is no reason whatsoever 
to take for granted that (|33|) defines a physical distance between two points 
of space, instead of, say, an optical distance. And if the distance is changed 
to a new one d(x\, x 2 ) which makes of B a point then the problem of choosing 
an appropriate radial coordinate has to be addressed concomitantly with the 
new space structure. 

We defend indeed in this paper the point of view that ( |3"3"D is an optical 
distance and that the physical distance between two points of space in the 
static frame of reference is instead 



X2 



d (xi, x 2 ) = / ds, (34) 
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where ds is the line-element of the principal transform of (^) . This leads of 
course to the conclusion that $ and \I/ in ( ^3|) are respectively the radial and 
tangential principal velocities of the speed of light 

To find the principal transform of the exterior Schwarzschild solution we 
shall start with its line-element written in curvature coordinates fl28|) . The 
general solution of Eq. ([H]) is then [13|, [14| 



R(r) = Qifi(r) + Q 2 f2(r), (35) 
Qi and Q% being two arbitrary constants of integration and 

/.W=r-^, / I (r) = jr^fr-?V (36) 
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Requiring 
leads to 



lim B! = 1 (37) 



Q 2 = l-Qi. (3* 



As we shall see in Sect. f|, the remaining constant Q± will be fixed by matching 
the exterior Schwarzschild solution to an interior one. Notice however, as it 
was already pointed out in the remarkable fact that if Qi = then the 
domain of staticity of the exterior Schwarzschild solution is in quo-harmonic 
coordinates the interval R G ]0, oof as for Schwarzschild's historical form. On 
the other hand quo-harmonic coordinates are intimately related to principal 
transformations and there is no doubt that in the sense of the 1-h-s of ( |2"2"| ) 
R = corresponds intrinsically to a point. 

The fate of the status of the Schwarzschild singularity remains therefore 
suspended to whatever we can learn from matching the exterior solution to 
an interior one using a global system of quo-harmonic coordinates. 



A new perfect fluid, spherically symmetric, 
static model 



We consider now the field equations (EB) under the general assumptions of 



staticity and spherical symmetry which led to the line-element ([[]), to which 
we shall add as a simplifying assumption the constancy of the energy density: 

p = constant. (39) 



A simplified and incomplete version of this point of view was developed in 1 11 and |jl 
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Using our notations and curvature coordinates Eqs. (2"5|) can be written as 
the following system of three first-order differential equations: 

2B- 3 B' + r-- 1 (l - £~ 2 ) = Krp, (40) 

2B-' 3 ((iBy = Kr, (41) 
2AA' + 2B- l B'A 2 = nra- 1 A 2 B 2 , (42) 

where 

a = — . (43) 
P + P 

Before writing down solutions of these equations let us remind that we 
want to obtain solutions that can be matched to the exterior Schwarzschild 
solution. This demands the existence of a value of r, say ri, such that 

i) on this 2-sphere, say TZ, the pressure vanishes, 

Pi = 0, (44) 



ii) and the functions A, B are continuous across 1Z. From (28) we have to 
require then, with obvious notations, 



/ 2m 

A^ = Af = ] Jl-—, A[B; = AtBt = l. (45) 

We shall discuss the continuity of the derivatives latter on. 
The general solution of Eq. fl40| ) is 

B~ 2 = 1 - - f r 2 pdr = 1 - — fr 3 - r$) , r > r , (46) 
r Jr 3r v ' 

where r is an arbitrary constant that we will assume non-negative. We shall 
refer from now on to r as the center of symmetry of the configuration. That 
r can be understood as being a point and therefore as a center will be fully 
justified in Sect. |[ 

One has to put r = to obtain the interior Schwarzschild solution. Since 
this particular solution satisfies the regularity conditions 

B = 1, B> = 0, (47) 

it has become customary to require the same regularity conditions from any 
other physically meaningful solution. This is unjustified for several reasons. 

i) It is assumed implicitly that the range of the radial variable must be 
r G [0,oo [.0 

3 This assumption and ro ^ leads to the singular solution of Volkoff |l5| and 
Wyman 111. 
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ii) An intuitive meaning of r is accepted before knowing whether a global 
C l model can be completed using curvature coordinates. We shall see that 
it can not. 

hi) More general regularity conditions can be accepted because they do 
not contradict any basic mathematical or physical principle. 

iv) It ignores that other interpretations for the line-element (0) are avail- 
able that do not require these regularity conditions. 

The values of B and its derivative at the center r = r are now 

B = 1, B' = ^npr . (48) 
A more general statement for the derivative B' is that 

B ' = ^ B3 ( 2r3 + r o)- (49) 

proving that B is a monotonously increasing function. 

From ( f45|) and (^) we obtain the mass or the density depending on which 
parameter we may wish to consider as given: 

1 / o ,\ 6m 
m= r p{r 1 -r ), p= -^—^ . (50) 

Taking into account (fil ) the solution of ( |41| ) is 

oB = -B 1 --k P rB 3 dr. (51) 

P 2 Jr 



The values of a and o' at the center are 

P 2 Jr 

Taking again into account (EH) the solution of (B3) is 



1 1 /" r i 1 
(T = -Bi —k / r5 3 dr, cr = -«r Q (1 - a p) . (52) 

P 2 Jrn 2 



A 2 = B~ 2 exp ( — k f 1 rB 2 a- 1 dr ) . (53) 



The values of A and A' at the center are 



A = exp (-7;^ r ^ 2(J 1 rfr ) > = ^ Kr oPo^o- (54) 

Notice that neither a', nor A' are zero at the origin r = r$. This is again 
a mild departure from conventional requirements of regularity that does not 
contradict any mathematical or physical principle we are aware of. 
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Let us prove that the pressure p is a monotonously decreasing function 
of r. In fact from (51]) and (|49|) , using (f43|)-(f44D it follows that at ri we have 
&i = p _1 and 

°'x = (r\ - rg) > 0, (55) 

because r\ > r . Since a is increasing at r 1; it will be monotonously increasing 
along the whole interval [to,?"!] unless there is a point r < r 2 < ri where 
= and cr 2 < which is impossible because if a' 2 = from ( pT| ) and ( pE9"| ) 
it follows that 



_3r| 
2r| + r 



a 2 = 7T^- s °"1 > ( 56 ) 



This completes the proof of our statement. We shall see in Sect. 5 that 
p (and, thus, the pressure at any other point) is finite except in the most 
extreme configurations one can think of. We shall also discuss there when 
the condition p < p is satisfied. 

The behaviour of A and a are very simply related. In fact, from (|51|) 



and fl53|) its follows that A/a = C with C = pJl — 2m/r 1; the value of this 



constant being derived from ( fj4[) and (4q). 
Let us write r in units of 2m: 

r = 2m a, a > 0. (57) 

Then from ([)t]) it follows that the mass m > is the solution of the equation 

6m + npa 3 (2m) 3 = npr\. (58) 

If < a < 1 then, for any given density, the viable values of r\ > r are 
bounded from above. This is a well-known result when a = 0, corresponding 
to the interior Schwarzschild solution, that follows immediately from (|46|) 
but it holds true also when r$ belongs to this more general interval. In fact, 



if we assume that r\ > 2m, as it is necessary from (|45|) to guarantee that 
Bf 2 > 0, then from @ it follows that 

— > (1 - a 3 )r 2 , (59) 
up 

which proves our assertion. 

On the contrary, if a > 1 all values of ri are viable whatever the density. 
In fact, the values of m > derived from (|58|) are such that 2m < r\ because 
otherwise we would have 

6m < Kp{\ - a 3 )(2m) 3 < 0, (60) 
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in contradiction with the positivity of m derived from (|58|). 

The preceding remarks prove that r = 2m selects a distinguished class 
of solutions that allow the consideration of models with arbitrary density, or 
mass, with arbitrary large values of r±. 

Consider the sequences of models with different values of a > 1 that 
can be constructed with a given mass m and decreasing values of r±. The 
sequence 

r = 2m (61) 

is the only one that allows to reach the limit where r at the exterior cov- 
ers the interval ]2m, oof of the maximum domain of staticity of the exterior 
Schwarzschild solution. We shall see in the following sections that this cor- 
responds to the most extreme models with point-like sources. And therefore 
this is the sequence that allows to construct models with arbitrary mass, or 
density, with arbitrary large or small size. 

Another striking difference between the two cases a = and a > comes 
from the comparison of the geometry described by the 3-dimensional met- 
ric (fj). Using the orthonormal co-basis 

l = dr, 6 2 = rd6 6 3 = r sin 6 dtp, (62) 

the non-zero components of the Ricci tensor are 

Rn = |§(2r 3 + r 3 ), R 22 = R 33 = ^(4r 3 - r 3 ). (63) 

It follows from these formulas the well-known result that for tq = the three 
non-zero components are identical, meaning that the metric (0) has constant 
curvature. This is not the case if tq > whatever this value might be. This 
means that our solution is essentially different from the Schwarzschild one. 

We postpone the discussion of the physical acceptability of these solutions 
as well as the discussion about the geometrical nature of their center. 

4 Matching interior solutions to the exterior 
one 

We have already implemented the continuity of the functions A and B on 
r = r\. It so happens that there remains no freedom to add new conditions 
on the derivatives when using curvature coordinates both in the interior and 
the exterior solutions. In fact, the derivative of B is inescapably discontin- 
uous, and the derivative of A is already continuous without the necessity of 
requiring it. This is of course the same that happens when trying to match 



11 



the exterior and interior Schwarzschild solutions but let us remind how this 
comes about. 

We shall use again, when necessary, super-indexes ± to refer to quanti- 
ties that belong to the exterior or interior solutions, and the sub-index 1 to 
indicate that the corresponding quantity has been evaluated on the border 
of the object. From ( f40| ) and (45) we have 



2 jyi 

2BfB? + ^r = Kpn, (64) 



and from (E28|) we have 



2m 



2B^Bf + — = 0, (65) 
rf 



and therefore inescapably 

2Bf 3 (B+' - Bf') = -npn 0. (66) 
Also from (|2|)-(|45|) it follows that 

2A 1 A±' + 2Bi 3 B{' = npn, (67) 



while from (EST) we have 



2A 1 Af + ^ = 0. (68) 



Subtracting (|67f) from ( |64l) and using the last result proves as stated that 

2A 1 (A? - V) = 0. (69) 

Because of (^B[) the model thus constructed is not of class C l across the 
border of the object. Most authors of papers and books when dealing with 
Schwarzschild's solutions feel satisfied with an identical situation. This is 
a mistake because it perverts the nice axiomatization of General Relativity 
proposed by Lichnerowicz [17| and introduces an ambiguity in the theory 



if one wants to have an interpretation for the metric (0), and the coeffi- 
cients (|23|) , as we have for instance for the function A. Besides, it is well 
knownQ that when ( fPiJ) holds then an interior static spherical solution can 
always be matched to the exterior Schwarzschild metric. Below we prove it 
for our particular interior solution.^ 

4 See for example |J| or @ . 

5 Another example can be seen in [M. 
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Let us consider the differential equation fl2~l"|) with B 
Eq. (|46|). The expansion of B around r = r is 



B given by 



B 



1 + -npr (r 



r ) + O 



r ) 



(70) 



and therefore we have to distinguish two cases. 

i) If ro = 0, i.e., when the metric is the interior Schwarzschild solution 
then r = ro = is a regular singular point and the solutions that are finite 
at this origin are all proportional to the following one:[] 



R 



2gr 2 



arcsm 



in(y/gr) — ^/qrdl — qr 



q " r p - 



(71) 



These finite solutions are such that 



R-(0) = R =0, R-'(0) = R \ (72) 

the derivative at the origin Rq' remaining arbitrary at this point. 

ii) If r = ro > then the origin is an ordinary point of the differen- 
tial equation. Therefore we can choose arbitrarily the values of R and its 
derivative at this point: 

R-(r ) = R ~, R-'(r )=R -'. (73) 

Each solution R(r) of the differential equation defines new coordinates (|20|) 
which are quo-harmonic for the interior solution. But we do not have, and 
we shall not use, their explicit analytic expression. 

These initial conditions can be chosen such that if R + is the solution of 
the same equation (|21~|) , with B = B + given by Eq. (|28|) , corresponding to 
initial conditions on the border r = r\, 

R+ = R^ R£' = R^\ (74) 

then the global model defined by the coordinate transformation R = R + (r) 
at the exterior and R = R~ (r) in the interior is of class C 1 across the border. 
In fact let the new expression of the metric (|J) be 

ds 2 = B 2 dR 2 + D 2 R 2 dn 2 , (75) 

where 

B r 

B _SR' ° = R- (76) 

6 This case has been considered by P. Teyssandier Q . 
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B, D and D' are continuous because R and R' have been chosen to be 
continuous. From the equation (ETC) that satisfy both R^ 1 we can derive that 



Rf" - i?r" - s^r 1 [bV - sr'j ^ = °- ( 7? ) 

From (|76|) we have 



Therefore, 



5' = fi'i?'- 1 - BR'- 2 R". (78) 

_ = B'R'~ l = R'~ 2 B' - R'~ 3 BR" (79) 
dR 



and 



dB^ dZ^ 2 , r ... , 



dR dR 

Using (|77|) in this equation we finally obtain 

dBf dBl 
~~d~R dR ~ 



R[- 2 Bf - B^' - R' 1 ~ 6 B 1 Rf - . (80) 



0. 



which proves that the global quo-harmonic model is of class C 1 . 

Once Ri and Ri' have been found integrating (|2lD with initial condi- 
tions (|73|) the constants Qi and Q 2 of ( p5]) can be found solving the system 
of linear equations 

Ri = QMn) + g 2 / 2 (ri), R-,' = QMn) + Q 2 f 2 {r l ). (82) 

Since the two functions f\ and f 2 are independent this system has always a 
solution. 

In case i) above the supplementary condition (^) determines the re- 
maining arbitrary constant Rq' , while in case ii) this condition gives only 
a relation between the two initial conditions (|73|) . To obtain a second rela- 
tion and make our model completely determine let us consider the principal 
transform of the interior metric of space. By definition we shall have (|22|) 
with R = R~ and appropriate functions B~ and C~ . If we accept, as 

we have done, that the structure of space in the static frame of reference of 
the source is given by the 1-h-s of this equation and we want to have a center 
for our configuration and no hole in the space we are forced to choose our 
function R with the initial condition 

R~ (r ) = i? " = 0. (83) 

Since the functions R(r) and r(R) are not explicitly known when ro > 0, 
neither is known explicitly the space-time metric (HD using the radial coordi- 
nate R. Therefore it is not possible using R to discuss directly the behavior 
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of the causal geodesies that go through the center of the configuration. But 
this can be done looking first at this problem using the original curvature 
coordinate r. Let us consider a causal geodesic reaching P in with r = r at 
some time t. Since the space trajectory is contained in a plane, say II, of the 
auxiliary Euclidean space with line-element dr 2 + r 2 dQ 2 , we can use polar 
coordinates such that II is the plane 6 = tt/2. Let (p be the value of the 
azimuth angle at p n ; and let r = v r and ip — v v be the two corresponding 
derivatives. Let us consider now the causal geodesic starting at time t from 
the point P out defined by r = r and ip = ip + it with derivatives given by 
r = —v r and ip = —v v . The space-time trajectory of this geodesic branch is 
contained also on IT and is the symmetric image of the incident branch with 
respect to the diameter joining P in and P ou t- We see therefore that the pro- 
cess of contraction of the sphere r = r to the point R = described by (|20|) 
joins the end-point of the incoming branch of the geodesic with the origin of 
the outgoing one. And in the process the matched geodesic is as smooth as 
they were the two branches. Notice also that, since starting from any point 
there are always a bunch of geodesies reaching r = r$ with different values of 
9 and ip, our construction of the center of the new models endows this center 
with a focalising property. 

The quo-harmonic class of coordinates derived from the radial coordinate 
R are natural coordinates in the sense that they have the following properties: 
i) they exhibit the center as a point in the sense of and ii) they make the 
space-time metric smooth across the border of the spherical source. Another 
system of natural coordinates sharing these two properties are the Gauss 
coordinates derived from the radial coordinate: 

R G = f r Bdr. (84) 

Jro 

Instead for large values of r the quo-harmonic R behaves as 

R = r- in + 0(l/r 2 ), (85) 



while the Gaussian Rq behaves as 



R = r — m 



, ( m 
1 + In — 

V2r 



0(l/r), (86) 



which makes the asymptotic behaviour of the metric to depart a little bit 
more from the Newtonian intuition. 

In the next section we shall discuss the numerical solutions and in the 
final section we shall discuss the physical relevance of them. 
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5 Numerical analysis of the model 



The numerical study of the model is greatly simplified by using the following 
dimensionless variables: 

x = , a — Xq = — , b = X\ = . (87) 

2m 2m 2m 

Now, using fl50|) the expression fl46[ ) reduces to 

B~ 2 = l-- X ^j. (88) 
x b 6 — a 6 



5.1 The pressure 

For the interior solution corresponding to x G [a, b], Eq. ( |o"T| ) reduces to 




which can be written in terms of elliptic functions, but is far more easily 
computed numerically. 

Since the pressure decreases from the center r = tq to the border r = ri, 
to make sure that it does not go to infinity it is enough to check that it does 
not diverge at r = vq, i.e., that /io = does not become zero (or negative). 
This can be easily checked by numerical quadrature. In Fig. [I] we see that, 
in the case a = 1 (r = 2m), fi only vanishes for b = a, which corresponds 
to the limit case in which the interior reduces to a point. The same happens 
for a > 1, but for < a < 1 the pressure becomes infinite at some interior 
point unless the matching radius r% = 2mb is greater than a given value ri m i n , 
which is displayed in Fig. |2|. In the case r$ = corresponding to the interior 
Schwarzschild solution, the well-known minimum value is r lmin = |m. 

One may also check numerically that the dominant energy condition p < p 
—i.e., fi > fi > = \ii{b) — is satisfied always if r > 2m and r\ > r . 
For < r < 2m, however, it is satisfied only for values of ri greater than 
a minimum value r lmin which happens to be higher than the r lmin discussed 
above, as displayed in Fig. [2[ In the case of the interior Schwarzschild solution 
(r = 0) one has fi min = |m. 

Thus, we conclude that 2m is the smallest value of r for which the 
physical conditions on the pressure do not put any limit on the matching 
radius r\ > r . 
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5.2 Quo- harmonic coordinates 



To compute the quo-harmonic coordinates (|20|) we have to solve the differen- 
tial equation (|2l|) , which in the dimensionless variables of this section reduces 
to 

2xHS" + (AH - 2x 3 - a 3 ) S' - 4 (o 3 - a 3 ) 5 = 0, (90) 

where 

a prime indicates derivative with respect to x and we have defined 

H = x(b 3 - a 3 ) B~ 2 = x(b 3 - a 3 ) - (x 3 - a 3 ) > 0. (92) 

In the case a = (corresponding to the interior Schwarzschild solution), 
the origin is regular singular, but the solution satisfying S(a) = 0, S'(a) = 1 
exists in [0,b] and can be explicitly written [TJ| as (|7T|), which reduces to 



3b 3 r 



S 2x* 



& 3/2 

arcsin (V 3/2 x) - xy/l - b- 3 x 2 ] . (93) 



In the remaining cases (a > 0), the differential equation is linear, the 
coefficients are continuous, and H does not vanish for x G [a, b] (we are 
excluding the limit case b = a), so that there exists a unique solution defined 
for a < x < b that satisfies the initial conditions S(a) = 0, S'(a) = (3, where 
(3 is a constant to be computed later. In fact, S = (3S, if S is the solution of 
Eq. ( |90|) satisfying 

5(a) = 0, S'(a) = 1. (94) 
Now the numerical method to compute R is straightforward. After se- 



lecting the parameter a, we compute S by solving ( |90|) with (|94). Then the 
auxiliary condition (R3) and the matching conditions ([52]) read 



= Qi0i(6) + (1 - Qi) <&(&), /?£'(&) = Q^(6) + (1 - Qi) <? 2 (&), 

(95) 

where functions (BBl) are written as 









H) 







oi (x) = = x — ^(x) = = l/l ( x ). (96) 

y v J 2m 4' y v ; 2m V x V 4/ K ' 

Since Eqs. ( |9~5"D are readily solved for (3 and Q±, R is given by 

g _ _ f PS(x), for a < x < 6, 

" 2m ~ 1 Qigi(x) + (I - Qi)g 2 (x), for x > 6. 1 J 
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Two particular cases, for a = 0, 1 and 6 = 2 (i.e., for ro = 0, 2m and r\ = 4m) 
are displayed in Fig. [| Notice the different definition ranges of R in curvature 
coordinates. 

In Fig. |] the coefficient Q\ is displayed for a = 1 and different values 
of b. One may see there that in the limit b — > a = f , i.e., when the interior 
collapses to the point R = in quo-harmonic coordinates and to the center 
ro = 2m in curvature coordinates, Q\ vanishes, so that the radial coordinate 
R associated to quo-harmonic coordinates is given everywhere by R = i^O"), 
as pointed in Sect. [3[ 



5.3 Metric coefficients and principal transform 

The metric coefficient B is given, in the dimensionless curvature coordinates, 
by (|8"8j ) for a < x < b and by B — (1 — 1/x) for x > b. Similarly, the 
other metric component in curvature coordinates is given by Eq. (0), which 
reduces to 

(9£ 



A 2 = 5" 2 exp f- J\b 2 ^ 1 dx^j , 



for a < x < b and is A = B -1 for the exterior x > b. Since we can compute 
numerically the quo-harmonic coordinates, they can be used to display the 
metric of different models in the same physical coordinates. 
The functions in Eqs. ( |2"BD now are written as 

• = §• » = i (99) 

and can be readily computed. 

For instance, in Fig. [5] we show the metric coefficients A and B for two 
models with a = 0, 1 (i.e., r = 0,2m). In both cases the matching of the 
interior and exterior metrics happens at the same location in quo-harmonic 
coordinates: at the spherical surface of radius S = 2 (R = 4m), which in 
curvature coordinates has radius x ~ 2.7134 when a = and x ~ 2.7916 if 
a = 1. The continuity of the metric coefficients and their first derivatives is 
apparent in the figure. The functions $ and ^ corresponding to the same 
special cases are plotted in Fig. |[ 

6 Concluding remarks 

Among all possible values of ro two of them are distinctly distinguished: 
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i) the value r = because it leads to the most well-behaved models, 
although only in a restrictive range of the density p and the radius of the 
interior configuration r\ or Ri, and 

ii) the value r = 2m because this leads to the closest sequence to 
Schwarzschild's one, allows unrestricted values of the parameters p and r\ 
or Ri and contains the emblematic most extreme configuration we can think 
of limit, namely that with a point-like source. 

Whether or not new considerations or difficulties will suggest selecting or 
excluding particular values of r remains to be seen. But if only one new 
value remains then be pretty sure that this value will be r = 2m. 

There exists a range of parameters p and r\ for which Schwarzschild's in- 
terior solution and ours coexist. Nevertheless the two solutions are quite dif- 
ferent, whatever these parameters might be. This is so because the geometry 
described by (|2|) is homogeneous and isotropic in one case and inhomogeneous 
and anisotropic in the other. 

We have assumed for simplicity that the energy- density was constant. It 
might be objected that this is not a realistic "equation of state" to describe 
the properties of our objects, and in particular the extreme ones, i.e., those 
without analog based on the interior Schwarzschild solution. But no equation 
of state will be realistic, and would be un-realistic to guess one, as long as 
new physics is not available for them. Not to mention that to assume the 
existence of an equation of state is also a simplifying assumption, as it is to 
use a perfect fluid description. The main properties of our models do not 
depend crucially, we believe, on these simplifying assumptions. They depend 
instead crucially on the identification of r = 2m with the origin and the 
center of of symmetry of the model and also on the re-interpretation of the 
metric @ based on the consideration of principal transformations. 

We do not claim that such objects will be found in nature, but we do claim 
that if they are found General Relativity is a good theory to understand its 
geometrical properties and gravitational field. We also claim that it is a 
sound attitude to look for them. Everybody knows where there is a good 
chance to discover them. 
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Figure 1: Function /i for a = 1 and small values of b. 
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Figure 5: Metric coefficients A and B when the matching radius is R\ = 
4m and r = 0, 2m. The same R coordinate associated to quo-harmonic 
coordinates is used in both cases. 
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Figure 6: Functions <£> and \& when the matching radius is R\ = Am and 
r = 0, 2m. The same R coordinate associated to quo-harmonic coordinates 
is used in both cases. 
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